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ABSTRACT 

It has recently been suggested that conduction-driven magnetohydrodynamic (MHD) instabilities may oper- 
ate at all radii within an intracluster medium (ICM), and profoundly affect the structure of a cluster's magnetic 
field. Where MHD instabilities dominate the dynamics of an ICM, they will re-orient magnetic field lines 
perpendicular to the temperature gradient inside a cooling core, or parallel to the temperature gradient out- 
side it. This characteristic structure of magnetic field could be probed by measurements of polarized radio 
emission from background sources. Motivated by this possibility we have constructed 3-d models of a mag- 
netized cooling core cluster and calculated Faraday rotation measure (RM) maps in the plane of the sky under 
realistic observing conditions. We compare a scenario in which magnetic field geometry is characterized by 
conduction driven MHD instabilities to that where it is determined by isotropic turbulent motions. We find that 
future high-sensitivity spectro-polarimetric measurements of RM, such as will be enabled by the Expanded 
Very large Array and Square Kilometer Array can distinguish between these two cases with plausible exposure 
times. Such observations will test the existence of conduction-driven MHD instabilities in dynamically relaxed 
cooling core clusters. More generally, our findings imply that observations of Faraday RM should be able to 
discern physical mechanisms that result in qualitatively different magnetic field topologies, without a priori 
knowledge about the nature of the processes. 
Subject headings: galaxies: clusters: general - instabilities - magnetic fields - MHD - plasmas - polarization 



1. EVITRODUCTION 

Understanding the role and evolution of magnetic fields 
in clusters of galaxies is of significant importance for many 
questions including the origin of cluster magnetic fields, the 
interaction of active galactic nuclei (AGN) with the intraclus- 
ter medium (ICM), and physical processes operating within 
the ICM plasma. 

One of the key techniques used to obtain information about 
the strength and structure of cluster magnetic fields is the anal- 
ysis of Faraday rotation from polarized radio sources located 
behind and within clusters. The Faraday effect rotates the po- 
larization plane of electromagnetic waves traveling through a 
magnetized plasma, by an amount proportional to the (square 
of the) wavelength, plasma density, and the strength of the 
field in the direction of propagation. The sources' intrinsic po- 
larization need not be known, as the effect can be observed as 
a characteristic wavelength-dependent rotation measure (RM) 
signature. Observations of a few nearby clusters have estab- 
lished the presence of magnetic fields with typical strengths of 
few/iG in non-cool core clusters and in excess of 10 /xG in the 
centers of cool core clusters (Carilli ^ Taylor|2002[|Govoni &| 
|Feretti 2004; Fe retti & Giovan nini 2008). Detailed high res- 
olution RM images of radio galaxies in merging and cooling 
core clusters indicate that the RM distribution is character- 
ized by patchy structures of a few kpc in size. Furthermore, 
the RM distribution appears in general in good agreement 
with the RM expected for magnetic fields geometry character- 
ized by tur bulent motions with a power-law power spectrum 
(|Enffin & V ogt 2003 , Mu rgia et al'j2004| 



Lamg et al.||2008f|Bonafede et al.|2010 



Guidetti et al. 2008 



for e.g.). This base 
of knowledge will be greatly enhanced by radio observatories 
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such as the Expanded Very Large Array (EVLA) and Square 
Kilometer Array^ (SKA), which will provide the sensitivity 
to study both background and embedded sources with suffi- 
cient density to map out the RM patterns in many and more 
distant clusters. This will allow to study the structure of mag- 
netic field in "quiescent" regions of clusters, away from radio 
galaxies, as well as the effect of the radio-mode feedback in 
their vicinity. 

The interpretation of these results can strongly benefit from 
advances in the theoretical understanding of the dynamics of 
weakly magnetized, low density plasmas. The ICM is a di- 
lute plasma, i.e., the gyro-frequency of both the electrons and 
the ions is much greater than the collision frequency. The 
principal effect of this is highly anisotropic thermal conduc- 
tion that fundamentally alters the classical Schwarzschild cri- 
terion for convection ( |Balbus| 2000| l. Under these circum- 
stances, the outer regions of clusters where temperature de- 
creases with radius may become susceptible to a conduction- 
driven magneto-thermal instability (MTI). The properties of 



this instability have been studied in MHD simulations (Par- 
|rish & Stone||2005| [20071 jParrish et al.||2008l l. [Parrish et al. 



(|2008') found that MTI can profoundly affect the temperature 
distribution in the outer regions of a cluster in only a few bil- 
lion years, by re-orienting the lines of magnetic field to be 
preferentially radial, resulting in a thermal conduction at a 
high fraction of the Spitzer conductivity. 

The temperature profile in the centers of cooling-core 
galaxy clusters makes them stable to MTI. However, Quataert] 
(|2008) found that cooling cores are characterized by a sis- 
ter heat-flux buoyancy instability (HBI) that arises in regions 
where temperature increases with radius. An initial investi- 



gation of nonlinear evolution of the HBI (Parrish & Quataert 
|2008| l indicates that its primary effect is to re-onent the lines 
of magnetic field perpendicular to the temperature gradient, 

^ http://science.nrao.edu/evla 
^ http://www.skatelescope.org 



and thus strongly suppress heat conduction. The properties of 
HBI have been investigat ed in glob al 3-d models of isolated 
cooling core clusters (Parrish et al. 2009; Bogdanovic et al. 



and scale them to our assumed cosmology). 



2009 1. They suggest that, once the magnetic field lines have 
been wrapped into spherical surfaces surrounding the core, 
the effective thermal conduction is suppressed to a small frac- 
tion of the Spitzer value, leading to insulation of the core from 
further conductive heating, and to a subsequent thermal col- 
lapse. 



Very recently, Balbus & Reynolds (20101 discovered an as- 
sociated pair of overstabilities that affect precisely those con- 
figurations that are stable to the well-established HBI and 
MTI. They predict that configurations which tend to result 
from the non-linear evolution of the HBI have g-modes that 
are driven overstable by radiative loses. On the other hand, 
configurations which tend to result from the non-linear evolu- 
tion of the MTI have g-modes that are driven overstable by the 
conductive heat flux. The effects of these overstabilities for 
the ICM plasma thermodynamics and the properties of mag- 
netic field are yet to be understood. We do not consider them 
in this work and instead focus on MTI and HBI instabilities. 

The combination of MTI and HBI instabilities should lead 
to a characteristic structure of magnetic fields in some cooling 
core clusters. If these conduction-driven instabilities domi- 
nate the dynamics of the ICM, the magnetic field lines would 
be preferentially oriented radially in the outer region and az- 
imuthally within the cooling core. This may not be the case 
in all clusters - it has recently been shown that MHD instabil- 
ities can be overwhelmed by even moderate levels of driven 
turbulence (Ruszkowski & Oh 2010, Parrish et al. 2010), such 
as might result from sub-cluster mergers, motions of cluster 
member galaxies, or various forms of AGN feedback. Which 
(if any) physical mechanism dominates in the ICM depends 
sensitively upon the magnitude and distribution of turbulence, 
which is currently only poorly understood. Motivated by this 
question and guided by the MHD simulations, we have con- 
structed 3-d models of clusters for two distinct cases and sim- 
ulated the RM maps that might be observed with high sensi- 
tivity radio polarization measurements. 

This paper is organized as follows. We describe the ingre- 
dients of our cluster models in § [2] and the main properties 
of derived RM maps in § l3] We discuss the importance of 
depolarization and diffuse emission in § HI then present the 
discussion with conclusions in §|5] 



2. CLUSTER MODELS 

2.L Intracluster medium and magnetic fields 

Any magnetized cluster acts as a Faraday screen for polar- 
ized sources located behind a cluster or in the cluster itself. 
Modeling the RM signature of such a cluster will require sev- 
eral main ingredients: the density and temperature distribu- 
tion of electrons in the ICM of the cluster, the 3-d structure 
of the magnetic field, and the density and fluxes of polarized 
background sources. We construct models in a Cartesian co- 
ordinate system {x,y^z) with a cubic spatial domain defined by 
X = zLL, y = ±L, z = ±L, where L = 300 kpc. For the electron 
density and temperature distributions we adopt analytic ap- 
proximations based on the XMM-Newton observations of the 
Perseus cluster that capture the radial behavior of these two 
parameters (we use expressions from jChurazov et al. 2003 
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where r = (x^+y^ + z^y'^ is in units of kiloparsecs and r^ = 
56.5 kpc, rti = 197.7 kpc, and r^. = 70.6 kpc. The central 
electron number density and temperature are ne{Q) = 4.3 x 
10-2 cm-3 and Te(0) = 3keV, respectively. 

We consider two models of the cluster's magnetic field. 
Our first model (model A) is motivated by the results of re- 
cent MHD simulations of instabilities in the ICM, where the 
magnetic field is preferentially azimuthal within the cooling 
core and radial outside of this region. Hence, in a spheri- 
cal polar coordinate system (r, 9, </>), where 9 = Q is aligned 
along the z-axis, the field structure within the cooling core, 
r <rc = 200 kpc, is described in terms of only 9 and </) com- 
ponents ( [Bogdanovic et al.|2009[ ): 

fie = 2Bo(l + sin(27rr/ri))sin6'cos(20), (3) 

B^ = 2Bq{\ + sin(27rr/r2)) sin(36') - (4) 

Bo(l + sin(27rr/ri))sin(20)sin(26'), 

Boif) is chosen so that the value of the plasma parameter, 
/3 = SnngkTg/BQ = 100, is constant everywhere within the 
cooling core region. This implies Bo(0) = 7.3 /iG at the clus- 
ter center, the field strength in the range 2.3 x 10-^ -36.2 //G, 
and mean strength over the computational volume of 3.3 /iG. 
ri = 7.5 kpc and r2 = 24 kpc are the coherence lengths defin- 
ing characteristic radial scales on which magnetic field vector 
changes direction. The magnetic field thus changes direction 
8 to 26 times across the cool core radius. In addition to these 
two scales, equations |3] and |4] also capture the field geome- 
try expected to arise as a consequence of the HBI: field lines 
wrapped onto the spherical surfaces within the cool core. This 
implies that in model A magnetic field reversal also occurs on 
a range of spatial scales associated with the spheres of differ- 
ent radii and up to the size of the cool core, 2r^ = 400 kpc. It 
follows that 400 kpc is the maximum scale for magnetic field 
fluctuations in model A. 

The magnetic field structure outside of the cooling core, r > 
re, is described by a radial component. 



Br- 



BrO 

r2 ' 



(5) 



where B^o is defined from the condition B^(rc) = {B^irc) + 
B\{rc)). The simple analytic form for the field structure given 
by equations ([3])-(|5]l satisfies the condition cZ/vB = every- 
where except at r = r^, where the field lines are transitioning 
from azimuthal to radial. In a reaUstic case, this transition 
would be more gradual and also divergence free. We never- 
theless expect that, except around r sa r^, our model should 
capture the salient properties of the magnetic field of a cluster 
affected by MHD instabilities across a wide range of radii. 

In our second model (model B), magnetic field lines are 
randomly tangled reflecting a different physical scenario in 
which the field geometry is set by the action of isotropic tur- 
bulence. While this theoretical hypothesis is most likely ideal- 
ized, it is physically motivated and well rooted in theoretical 
practice, so we use it as a control case to model A. Follow- 
ing Roettiger et al. ( 1999 1, we initialize the field geometry by 
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TABLE 1 
Summary of model properties 



Model Observatory Scenario 



Ns ds I]|a,,.o| J2\"n.l\ Y^l^iial I]|Q>i,3| 

(kpc) J]kn.m| Z]l"n.m| Z]l"'i"i| Z]l"n."i| 





SKA 


Ih, low 


449 


28 


8.1% 


39.1% 


15.3% 


2.1% 




SKA 


Ih, high 


7117 


7 


17.3% 


45.1% 


15.4% 


1.9% 


A 


SKA 


lOOh, low 


2250 


13 


7.8% 


43.5% 


17.6% 


2.5% 




SKA 


lOOh, high 


89592 


2 


14.3% 


50.0% 


17.7% 


0.7% 




EVLA 


9h, low 


65 


75 


11.1% 


28.2% 


10.5% 


3.9% 




EVLA 


9h, high 


337 


33 


9.0% 


39.5% 


15.8% 


2.4% 




SKA 


Ih, low 


449 


28 


25.3% 


15.6% 


11.6% 


5.4% 




SKA 


Ih. high 


7117 


7 


11.9% 


12.1% 


9.4% 


5.6% 


B 


SKA 


lOOh, low 


2250 


13 


19.3% 


16.0% 


10.4% 


5.7% 




SKA 


lOOh, high 


89592 


2 


9.2% 


7.5% 


9.1% 


5.7% 




EVLA 


9h, low 


65 


75 


31.1% 


16.6% 


13.2% 


4.4% 




EVLA 


9h, high 


337 


33 


20.5% 


19.6% 


11.0% 


7.6% 



Model - A (B) con'esponds to the instability-dominated (turbulence-dominated) cluster model; Obser- 
vatory - one of the two radio-observatoiies considered in the paper; Scenario - see the text; W, - number 
of background polarized sources within 600 kpc X 600 kpc area; rf, - mean separation of polarized back- 
ground sources; ratios - fraction of the total RM intensity in circularly symmetric shapelet multipoles. 



defining a magnetic field potential in Fourier space of A(A;) = 
AqA:"", where each Cartesian factor Aq has an amplitude 
drawn from a Gaussian distribution and a random phase, as- 
suring uncorrelated modes. We convert this into real space via 
a three-dimensional fast Fourier transform (FFT). The tangled 
magnetic field is then calculated as B = V x A. We adopted 
a = 17/6, which results in the Kolmogorov-like power spec- 
trum B^ ex fe^^'~"' ex k~^^'^. The smallest and largest magnetic 
structures produced by this power spectrum have the scales of 
Amin = 27r/^,nax = 7.5 kpc and A^ax = 27r/^min = 600 kpc, re- 
spectively. The magnitude of magnetic field is normalized in 
such way that the azimuthally averaged magnetic energy com- 
prises 1% of the thermal energy of the gas at all radii (i.e., 
[3 = 100, same as in model A). This implies the mean mag- 
netic field strength decreasing with the distance from the clus- 
ter center and Bo(0) = 7.3 /iG. The field strength in model B 
varies between 8.3 x 10"**- 29.8 /iG and its mean amplitude 
over the computational volume is 1 .3 /iG. 

2.2. Observational scenarios 

We place the model cluster at the redshift of the Perseus 
cluster, which is a suitable prototype for a nearby cooling core 
cluster. A t the distance of NGC 1275 (z = 0.0176; [Strauss et| 
al. 19921, the central bright galaxy in Perseus, 1" corresponds 
to 353 pc.^ To simulate future, high sensitivity RM surveys, 
we adopt the planned capabilities of the soon to be fully op- 
erational EVLA and the next-generation radio interferometer, 
5XA, at 1.4 GHz. 

EVLA consists of 27 25-meter diameter antennas that will 
provide an order of magnitude increase in sensitivity above 
the existing Very Large Array after it is upgraded with new 
receivers and electronics. Together with improved resolution 
and imaging, these new capabilities make EVLA an important 
tool for studies of Faraday rotation in the near future.^ Its 
technical specifications include a continuous coverage with 
full polarization capabilities between frequencies of 1 and 
50 GHz, a field of view (FOV) of - 0.25 sq. deg at 1 .4 GHz, 
and angular resolution as high as 6'evla = 1-3" (achievable 
in the largest array configuration; A-configuration) also at 

' Throughout this paper we adopted a ACDM cosmology with Hq = 
71 kms-' Mpc-' , n,„ = 0.27, and Oa = 0.73. 

' Currently, EVLA is on schedule for completion at the end of 2012. 



1.4 GHz (A = 21 cm). Furthermore, observations of the con- 
tinuum emission with the EVLA are expected to achieve the 
rm.s. noise limit of 1.6/iJy/beam between 1 and 2 GHz in 
about 9h of exposure. In the next step we will use this sensi- 
tivity limit to estimate the density of the background polarized 
sources on the sky that will be seen by the EVLA at 1 .4 GHz.^ 
With further improvements in sensitivity and survey speed 
over current instruments, SKA will be ideally equipped to 
study the origin and evolution of cosmic magnetism in the fu- 
ture dCarilli & RawUngs|2004l [Dewdney et al.|200"9HKrause| 



et al. 



2009| l.^ Planned technical specifications for the SKA 



include a square kilometer collecting area, continuous fre- 
quency coverage from 70 MHz to 25 GHz, a FOV of 1 sq. deg 
at 1.4 GHz, and angular resolution better than ^ska = 1" at 
the same frequency ( Schilizzi et al.|2007t[Taylor|2008| ). Note 
however that the frequency range for RM studies with the SKA 
will be relatively wide (^ 0.3- 10 GHz), and that the exact 
FOV and angular resolution depend on the frequency. Obser- 
vations with the SKA are expected to achieve an r.m.s. noise of 
0. 1 //Jy/beam area and 0.01 /iJy/beam area at 1 .4 GHz within 



Ih and lOOh of integration, respectively (CarilU & Rawlings 
[2004l l. 

Given the EVLA and SKA sensitivity limits, we follow ap- 
proach in Stepanov et al. ( 2008 ) to estimate the number den- 
sity of polarized bac kground sources at 1.4 GHz expected 
for our model cluster. [Stepanov et al. ( 2008 1 extrapolate the 
source counts at 1 .4 GHz from P = Q.5 mJy to the limiting flux 
density of Pmin = 0.05 fiJy using a power-law relation between 
the cumulative number of sources per square degree and po- 
larized flux density, Nn{> P) oc P~'>'. The exponent 7 can take 
a range of values between 0.7 and 1.1, which correspond to 
the low and high bounds for the source counts at a given flux 
limit at 1 .4 GHz. We consider exposure times of 9h in case 
of ^yZA and Ih and lOOh for SKA, both in the low and high 
source count scenarios, thus obtaining a total of six different 
observational scenarios for both models A and B. In Table [1] 
we list the total number of background polarized sources (A^^) 
with polarized flux densities > Sap above the noise, within a 
600 kpc X 600 kpc area centered on the cluster. Also shown is 

** Information about the EVLA capabilities and specifications was obtained 
from http://www.vla.nrao. edu/astro/guides/vlas/cun'ent/ . 

' At this point a commencement of Tiill science operations with the SKA is 
planned for 2020. 




Fig. 1 . — Simulated continuous Faraday rotation measure maps for an instability-dominated (model A) and turbulence-dominated (model B) cluster. Panel size 
is 600 kpc on a side and the color scale shows log |RM|. In model A, two distinct regions are visible: the magnetic field is dominated by HBI inside the cooling 
core, and MTI outside. The diagonal feature at large radii is an artifact of the model. In model B, the cluster's magnetic field is randomly tangled, producing a 
patchy RM distribution. 




Fig. 2. — Faraday rotation measure maps simulated for a 9 hour exposure with EVLA in a high source count observational scenario. Panels show instability- 
dominated cluster considered in model A (left) and turbulence-dominated cluster from model B (right). Panels show the same region and color scaling as in 
FigurefT] 



the mean separation between the sources (c/j) at the redshift 
of the cluster 

Note that in the estimation of the number of background po- 
larized sources we simply extrapolated the number of sources 
from 0.5 mJy down to the microJy and nanoJy flux densi- 
ties. The nature of the dominant population of sources that 
EVLA and SKA will detect at such low flux densities is un- 
certain, as they lie beneath the surface of the deepest radio 
surveys that are currently available. It has been proposed that 
microJy and nanoJy sky i s dominated by radio-quiet AGN 
and star- forming galaxies (Jarvis & Rawlingsl|2004 Wilman 



et al. 2008)1 as well as faint ellipticals and dwarf galaxies 



(jPadovani 1120101). It is also worth pointing out a simplifica- 
tion in our modeling which stems from an assumption that all 
background polarized sources are point like, while in reality 
some low redshift background sources will have a resolved 
extended structure. We nevertheless expect these to comprise 
a small fraction, as a majority of sources are expected to re- 
side at high redshifts. 



3. FARADAY ROTATION MEASURE MAPS 
3.1. Construction and properties 

The polarization of electromagnetic waves traveling 
through the cluster's magnetic field is rotated by an angle 
Ax = RM)? /(\ +z)^, where A is the observing wavelength 
and z is the redshift of the Faraday screen. We calculate the 
effect of this Faraday rotation for our model cluster by numer- 
ically integrating the expression RM= 812radm~^ Jn^B dl, 
where n^ is the electron number density in units of cm~^, B 
is the vector of magnetic field in /iG, and 1 is the vector de- 
marking the depth of the magnetic screen in kiloparsecs as 
measured along the line of sight of an observer placed at in- 
finity, at an arbitrary orientation with respect to the cluster. 
The result is a two dimensional, continuous map of RM in- 
tensity with ^1" resolution shown in FigurefT] We convolve 
this map with a map of point like, randomly distributed po- 
larized background sources, where the number of sources in 
the map is Ns- For visual representation in Figures|2]|3] and|4] 
the convolved maps are smoothed with a box function, where 
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Fig. 3. — Simulated SKA Faraday rotation measure maps for an instability-dominated (model A) cluster The four observing scenarios include short or long 
exposure times, with low or high background source densities. Panel size and color scale are the same as in previous figures. 



the size of each box is chosen so that it contains 5 polarized 
background sources on average. This implies smoothing box 
of the size h^ox = 69 kpc for the EVLA 9h, high scenario and 
hhox = 60, 15, 26.4, and 4.2 kpc for the SKA scenarios Ih, low, 
Ih, high, lOOh, low, and lOOh, high, respectively. We list the 
properties of the RM maps in Table [1] 

Figure l2] shows maps for the models A (left panel) and B 
(right panel) calculated for a 9h exposure with EVLA in the 
scenario with an assumed high count of background polar- 
ized sources. The RM structure in both appears marginally re- 
solved. The 9h, low scenario results in unresolved RM struc- 
ture (the RM map for this case is not shown but its properties 
are listed in Table [1] for comparison). This is not surprising 
given the low number density of background sources within 
the considered area with average spacing between the sources 
of ds ~ 75 kpc, demarking the effective size of the "resolution 
element" in this scenario. Figures[3]and|4]show the RM maps 
for models A and B, realized in each of the four observational 
scenarios with the SKA. In both, the instability-dominated and 
turbulence-dominated scenarios, the quality and spatial sam- 
pling of the maps increase with the density of background 
polarized sources, to the degree that the most optimistic ob- 
servational scenario with the SKA, lOOh, high, almost exactly 
replicates the features of the continuous RM maps shown in 
Figure [T] This is implied by the hierarchy of characteristic 



scales, 6'sKA < ds < Amin, meaning that both the angular res- 
olution of the SKA and the density of the background polar- 
ized sources are sufficient to fully capture the spectrum of the 
~kpc scale RM variations in this model. 

In all observational scenarios of model A it is possible to 
discern two distinct regions, which are characterized by the 
HBI (core) and MTI (outer region) instabilities. The dark 
swirl patterns in the cooling core region, with low values of 
RM, arise where the magnetic field changes direction along 
the line of sight. Because of the multiple magnetic field re- 
versals along the lines of sight, the RM intensity ridges in 
the core do not correspond to the magnetic field structure in 
a trivial way. Nevertheless, they carry an imprint of the az- 
imuthal distribution of the field lines, and if indeed present in 
cooling core clusters, may be one of the characteristic features 
to search for in RM surveys of clusters. Outside of the core, 
where the magnetic field geometry is radial with decreasing 
magnitude, the RM intensity also decreases smoothly with ra- 
dius. The diagonal feature apparent in maps calculated for 
model A, in the left panel of Figure |2] and all panels of Fig- 
ure l3] is an artifact of our model: the underlying symmetry in 
the radial component of the magnetic field causes a cancella- 
tion in the RM along these lines of sight resulting in low val- 
ues of RM. The black specs occur in places where the number 
of background sources per smoothing box falls to zero. The 




Fig. 4. — Simulated SKA Faraday rotation measure maps for a turbulence-dominated (model B) cluster, in four observational scenarios. Panel size and color 
scale are the same as in previous figures. 



lower cutoff value applied in maps (lOradm" ) is compara- 
ble to the maximum error of RM reached at t he limiting flux 
densit y f min = 5ap with the EVLA and SKA ( Stepanov et al. 
[2008] l. 

In model B, shown in the right panel of Figure [2] and Fig- 
ure |4j the RM patterns exhibit a noticeably different, patchy 
distribution. The local maxima (minima) of the RM distribu- 
tion are again associated with regions where magnetic field 
reversal results in an enhancement (a cancellation) along the 
line of sight. The main difference in Model B is that the can- 
cellation effect is more pronounced relative to model A, due 
to the tangled geometry of the magnetic field, even though the 
assum ed m agnetic field strength is comparable in both models 
(see § |2.1| i. This results in the maximum RM value about an 
order of magnitude lower than in the model A, where the mag- 
netic field lines exhibit uniformity and azimuthal structure on 
large scales. 

3.2. Discrimination between cluster models 

An image analysis technique well suited to characterizing 
the different geometric patterns in the RM maps is polar 
shapelets (Massey & Refregier 2005). This is based on the 
unique decomposition of localized objects into a series of or- 
thogonal basis functions that explicitly separate modes with 
different rotational symmetries. Useful forms exist in both 2-d 



and 3-d. The former can generally be used for accurate object 
photometry and astrometry ( |Kuijken|[2008| l, as well as mor- 
phological classification of the im ages of galaxies ( [Kelly &| 
McKay 2005''Masse y et al.|200 7), magnetograms of sunspots 
1 Young et al. 2005), and the response of the human visual 
cortex ( [Victo r et al. 2009). Its convenient mathematical prop- 
erties and intuitive interpretation also make it a particularly 
effective morphology estimator for clusters' magnetic fields. 
We fit each image /(x) in Figures [2] [3] and |4] as a weighted 
sum of shapelet basis functions Xn mix^ sucnthat 



30 n 
n=0 m=-n 



iXn,m{x,y), 



(6) 



where the (complex) coefficients fl„ ,„ describe the power in 
modes with n radial oscillations and m-fold rotational symme- 
try, similar to a localized Fourier transform. Higher «-orders 
also capture structure at increasing distances from the cluster 
core. In general, the sum over n can extend to infinity, al- 
though in practice we have arbitrarily truncated it to 30. The 
sum over m need include only every other term, because the 
intervening basis functions are explicitly zero. We focus our 
analysis on the inner 600 kpc region of the modeled cluster, 
which encloses the most interesting RM patterns, and show 
the magnitudes of the derived shapelet coefficients a„ ,„ in Fig- 
ures [51 |6l [71 
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Fig. 5. — Polar shapelet decomposition of the RM images for model A (left) and model B (right) in observational scenario 9h, high with EVLA. The amplitude 
of multipole coefficients is indicated by the linear greyscale, with darker colors corresponding to higher amplitudes. The (complex) shapelet coefficients also 
have phases, which indicate the orientation of each multipole, but these are not shown. 



The analysis of the maps for different observational scenar- 
ios within model A (or model B) give qualitatively similar re- 
sults, but with varying ratios of signal to noise. The key idea 
is that the shapelet decomposition efficiently isolates modes 
of rotational symmetry present in the RM maps, even in the 
presence of significant observational noise. Indeed, the dis- 
tribution of power in the two models is strikingly different. 
The shapelet decomposition of scenarios in model A shows 
almost all power within the range m = {-3,3}, while it fans 
out uniformly over m in turbulent model B. The radially sym- 
metric (m = 0) modes describe the average, uniform level of 
the RM intensity. The restriction of model A's deviations from 
dipole {\m\ = 1), quadrupole (\m\ = 2) and sextupole (\m\ = 3) 
modes with low orders of rotational symmetry can be under- 
stood in light of the equivalent 6 and cj) terms in equations (|3]l 
and (|4]). For each model, Tableflllists the fraction of power in 
shapelet coefficients with \m\ = to 3, as estimators similar to 



the shapelet asymmetry estimator (equation 61 of Massey et 
|al.|2007 D. 

In all the observational scenarios that we have considered, 
the magnetic field geometry of model A is distinguished from 
that of model B as a large fraction of RM power in modes 
\m\ < 3 and a decreasing power in higher-|m| modes. In sce- 
narios with a low density of background polarized sources, 
namely 9h, high (left panel of Figure l5]l and Ih, low (top left 
panel of Figure |6|, the observed RM measure map has in- 
sufficient resolution to contain all the \m\ = 3 modes, so the 
cutoff is present but less distinct. That the two EVLA and 
SKA observational scenarios are indeed similar is indicated by 
the comparable polarized source densities and similar distri- 
butions of the RM intensity (TablefT). Most excitingly, even in 
observing scenarios with noisier data, which may be achieved 
in longer exposures with EVLA and relatively short expo- 
sures with SKA, the large-scale magnetic field patterns are 
efficiently captured by polar shapelets. In the remaining ob- 
servational scenarios for the SKA, characterized by the higher 
density of background polarized sources {Ih, high, lOOh, low 
and lOOh, high), the distribution of power remains very ro- 
bust, as illustrated by the polar shapelets decomposition. 

In model B, the power is more uniformly distributed among 
the azimuthal modes (m-modes) while the distribution across 



the radial modes (n-modes) changes with the number density 
of background polarized sources. In observational scenarios 
with the lower density of sources (9h, high, Ih, low and lOOh, 
low), a significant fraction of the power is in lower n-modes 
and thus, the RM structure is captured on larger scales and 
not captured on small scales. With the increasing density of 
polarized sources some of the power shifts uniformly towards 
higher m and n-modes (see Ih, high scenario), indicating that 
finer RM structures are beginning to be resolved in this sce- 
nario. 

A practical implementation of the shapelet decomposition 
also requires advance estimates of a cluster's center and size, 
so that the basis functions can be constructed at a given loca- 
tion. For this analysis, we iteratively optimized the size of the 
basis functions to minimize residuals between the simulated 
RM map and its shapelet model. The scale factor of the basis 
functions represents the size of the Gaussian used in image de- 
construction. In combination with the maximum (truncation) 
value of n in the model, it determines its resolution in such 
way that smaller values of the scale factor allow higher res- 
olution models. The same scale size also sets the maximum 
spatial extent of the model, and there is a balance between 
the ability to model the large scale features and the high fre- 
quency detail, a choice that is optimized via x^ minimization. 
The used sizes of the basis functions thus naturally arise to be 
different in the various observational scenarios and between 
models A and B, depending on the detail of the RM patterns. 
This approach allows to constrain the characteristic scale of 
the smallest (resolved) RM patterns using this iterative proce- 
dure, however, our results are robust to setting a fixed physical 
scale size to the reconstruction. Also for the purposes of this 
analysis, we assumed that the center of a cluster would be 
known a priori. It would be possible to determine the best-fit 
center via iteration on the RM image itself, but it will likely 
be known to better accuracy in practice from independent (e.g. 
optical or X-ray) observations. Spurious offsets in the center 
primarily shift power in the dipole |m| = 1 shapelet modes. 
Indeed, small-scale structure near the cluster core in this real- 
ization, amplified by the high electron density, is responsible 
for the large ai,±i and 03, ±3 coefficients, which correspond to 
basis functions with a small spatial extent. If the center were 
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Fig. 6. — Polar shapelet decomposition of the RM images in model A (instability-dominated cluster) for four different observational scenarios with the SKA. 
Color scale is the same as in Figurels] 



determined from this data alone, this structure would pull the 
center around, s hifting some pow er between adjacent shapelet 
coefficients (M assey et al.|2007 1. 

We emphasize that the exact distribution of power in de- 
composed RM images is model-dependent and that neither of 
our two models should be regarded as a strict prediction of fu- 
ture observations. More generally, our analysis demonstrates 
that physical mechanisms that qualitatively adjust the mag- 
netic field distribution in clusters can be easily distinguished 
using observed RM maps, with plausible exposure times. 

4. DEPOLARIZATION AND DIFFUSE FOREGROUND EMISSION 

While the sensitivity of the instrument and (consequently) 
the density of observable polarized background sources 
strongly affect the precision of Faraday rotation measure ob- 
servations, the ability to tie the observed properties of the RM 
maps to the underlying magnetic field also depends on the ef- 
fect of depolarization. Depolarization is a reduction of the ob- 
served degree of polarization which may arise within the clus- 
ter itself due to its spatial extent (internal depolarization) or 
due to limitations in instrument capabilities (beamwidth and 
bandwidth depolarization). We now disucss the importance 
of beamwidth and bandwidth depolarization in the context of 
our models. 



Beamwidth depolarization arises when the minimum mag- 
netic field coherence length is smaller than the beam size 
of the radio instrument, and cancellation of the RM occurs 
within the beam. Given that the assumed beamwidth in our 
calculation is 1" = 353 pc and the magnetic field coherence 
lengths are much larger (> 7.5 kpc), no beamwidth depolar- 
ization is expected to occur in these scenarios. A similar ef- 
fect that can lead to a small loss of information is smoothing 
that we apply to all of our RM images. In our simulated ob- 
servations, the smoothing length varies between ^ 4 kpc and 
^ 70 kpc, depending upon the density of background sources 
in a particular observational scenario, and this determines the 
size of the smallest RM patterns that can be inferred with 
confidence from a smoothed image. However, this does not 
represent a fundamental limitation, since more sophisticated 
smoothing methods and scales could be adopted on real data. 

Bandwidth depolarization is the cancellation and averag- 
ing of the Faraday rotation measure that arises in polarimeters 
operating in fixed, wide frequency bands. Bandwidth depolar- 
ization will be largely eliminated by the EVLA and SKA, by 
acquiring RM measurements via the rotation measure synthe- 
sis method. RM-synthesis is based on multichannel spectro- 
polarimetry to enable the detection of weak, polarized emis- 
sion (Brentjens & de Bruyn^2005b). Most importantly, this 
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Fig. 7. — Polar shapelet decomposition of the RM images for model B (turbulence-dominated cluster) for four different observational scenarios with the SKA. 
Color scale is the same as in previous figures. 



also allows the simultaneous observation of a range of differ- 
ent RM values, and the separation of RM components from 
distinct regions (such as foreground and background struc- 
tures) along the line of sight. Indeed, the RM signal due 
to Faraday rotation of polarized light from real background 
sources is complicated by diffuse, intrinsic polarized emission 
from both the ICM plasma itself, and our own Galaxy. 

Diffuse, polarized foreground emission from cooling core 
clusters can be attributed to mini-halos, steep spectrum ra- 
dio sources associated with the ICM around a powerful cen- 
tral radio galaxy. The prototypical example is a ^ 450 kpc 
mini-halo at the center of Perseus. Brent jens] (2011 1 find 
that the diffuse emission observed at 350 MHz in the direc- 
tion of Perseus seems not to be related to the mini-halo but 
rather to the foreground emission from the Milky Way. Our 
Galaxy interferes with measurements by Faraday-rotating any 
extragalactic polarized signal and by adding its own polarized 
emission. RM of Galactic origin is typically ^ lOradm"^ but 
can be as high as 300radm~^ for objects close to the Galactic 
plane fSimard-Normandin et al.|1981|l. In the direction of the 
Perseus cluster for example, 'Brentj ens| ( |201 \\ measure a rela- 
tively low RM contribution (compared to our modeled values) 
in the range -50 to +100radm~^, which can in principle be 
disentangled from the RM map of the cluster using the RM- 



synthesis technique and is not expected to produce an RM sig- 
nal competing with that from the cluster. Thus, we conclude 
that the RM features modeled in this work are not expected to 
be significantly affected by depolarization and should domi- 
nate in magnitude over any component of RM contributed by 
the sources of diffuse polarized emission. The primary fac- 
tor that determines the efficacy of the RM maps in probing 
the magnetic field structure is the sensitivity of the spectro- 
polarimetric measurements. 

5. DISCUSSION AND CONCLUSIONS 

Even if MTI and HBI instabilities operate uninhibited in 
real clusters, the radial-azimuthal field geometry is likely to be 
perturbed by intermittent phases of AGN activity and merg- 
ers. As these perturbations are attenuated over time, they 
can source a new cycle of MHD instabilities. The relics of 
such events, including shocks, bubbles, and ridges of mag- 
netic field lines swept by intracluster galaxies, may also be 
recognizable in high quality RM maps - testifying to the past 
evolution of a cluster But with a characteristic time scale 
for the saturation of MHD instabilities of only a few billion 
years, the magnetic field lines can be driven to recover their 
orientation relatively quickly. As a consequence, a fraction of 
clusters may exhibit preferential magnetic field geometries. 
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despite episodic disruptions. 

Circumstantial evidence in support of this hypothesis may 
already exist in several forms. Recent detailed observations 
of the Virgo cluster revealed an isothermal region with re- 
markable azimuthal symmetry cocooned in the cool core of 
the cluster between the active radio-lobes of M87 (Mi llion et| 
5T]|2010|l. The metallicity of the isothermal gas in the same 



line of evidence comes from the study of magnetic field struc- 
ture ar ound galaxies in the Virgo cluster. |Pfrommer & DursT] 
(20101 find that in Virgo, which seems to be in transition to a 



region is non-uniform and clumpy. It is reasonable to expect 
that the high metallicity parcels of gas were uplifted from the 
low entropy, cool core region of the cluster, where the den- 
sity of stars is highest. That the parcels now have the entropy 
and temperature of the surrounding ambient plasma, but retain 
clumpy metallicity distribution, suggests that turbulent mix- 
ing was not efficient here and that the gas was instead heated 
by conduction. The presence and shape of the isothermal front 
indicate that the magnetic field geometry in this region may 
be predominantly azimuthal. Furthermore, the absence of ef- 
ficient turbulence indicates that HBl may operate unhindered 
in this region and that the instability may be more robust than 
suggested by some recent theoretical works ( [Ruszkowski &| 
|Oh 2010; Parri sh et al.||20 1(3) 1. However, at least one more 
physical phenomenon may produce a similar magnetic field 
topology: [Ruszkowski & Oh] ( |2010[ ) find that weak turbulent 
motions lead to trapped g-modes and result in gas motions 
that are preferentially tangential. They suggest that magnetic 
fields in clusters experiencing such g-modes can in principle 
become tangential even in the absence of thermal conduction 
and the HBl. While some thermal conduction seems to be im- 
plied in the case of the Virgo cluster, the concurrent presence 
of g-m odes ca nnot be eliminated. 

Juet t et al.| (J2010) draw attention to what may be more 
circumstantial evidence for MHD instabilities. As much as 
20% of the sample of 70 clusters presented by |Snowden et| 
[aL] ( |2008| l have puzzling temperature profiles that appear to 
be quasi-isothermal at ^Mpc radii ^''. The clusters show no 
signs of interaction, which could possibly have accounted for 
their temperature distribution, and they appear dynamically 
relaxed. Moreover, it is unclear whether models based on the- 
oretical studies of cluster properties could account for this un- 
usual class of objects (Nagai et al. 2007 1. If the class is shown 
to be unexpected, it may point to some missing physics in 
cosmological simulations. A possible explanation for the tem- 
perature structure of these objects is that thermal conduction 
operates very efficiently in their outer regions, an effect that 
arises as a natural consequence of the MTI instability. Ra- 
dial configuration of magnetic field lines in outer regions of 
such clusters could only persist if they are isolated and unper- 
turbed for sufficiently long periods of time, consistent with the 
observed properties of the quasi-isothermal clusters. Another 

'" The classification of the sample of quasi-isothermal clusters according 
to the thermal state of their cores has not been reported thus far. 



cool core, the global magnetic field has a predominantly radial 
orientation at large radii, which again suggests the operation 
of MTI in its ICM. 

In this study, we have evaluated the effects of two different 
physical mechanisms on the Faraday rotation measure of a 
magnetized cooling core cluster, in the context of the planned 
capabilities of the EVLA and SKA radio observatories. We 
compare a theoretical scenario in which conduction-driven 
MHD instabilities dominate the dynamics of the ICM, to a 
scenario in which magnetic field topology is defined by tur- 
bulent motions. We employ the polar shapelets image anal- 
ysis method to efficiently detect patterns in the RM image 
with specific rotational symmetries, and thus classify their 
morphologies. Within the bounds of our simple models we 
find that the two mechanisms can produce strikingly differ- 
ent RM patterns and that future spectro-polarimetric measure- 
ments will have sufficient sensitivity to discriminate between 
them that can be achieved in longer exposures with EVLA and 
relatively short exposures with SKA. We propose that the ef- 
fect of the HBl and MTI instabilities be sought for in dynam- 
ically relaxed cooling core clusters, and especially in the sub- 
class of clusters with quasi-isothermal temperature profiles 
at large radii. More generally, it should be possible to dis- 
cern physical mechanisms that result in qualitatively differ- 
ent magnetic field topologies from observed Faraday rotation 
measure maps, without a priori knowledge about the nature 
of the processes. Such observations will enable detailed in- 
vestigations into the behavior of MHD instabilities and other 
associated physical phenomena, which are of far reaching im- 
portance to a number of fundamental questions related to en- 
ergy transport in clusters. 



We thank the anonymous referee for thoughtful comments 
which helped to significantly improve this manuscript. T.B. 
would like to thank Tracy Clarke and Richard Mushotzky for 
stimulating discussions and useful suggestions. Support for 
T.B. was provided by the National Aeronautics and Space Ad- 
ministration through Einstein Postdoctoral Fellowship Award 
Number PF9-00061 issued by the Chandra X-ray Observatory 
Center, which is operated by the Smithsonian Astrophysical 
Observatory for and on behalf of the National Aeronautics 
Space Administration under contract NAS8-03060. T.B. and 
C.S.R. acknowledge support from the NSF under grant AST- 
0908212. R.M. is supported by STFC Advanced Fellowship 
PP/E006450/1 and FP7 grant MlRG-CT-208994. 



REFERENCES 



Balbus, S. A. 2000, ApJ, 534, 420 

Balbus, S. A., & Reynolds, C. S. 2010, ApJ, 720, L97 

Bogdanovic, T., Reynolds, C. S., Balbus, S. A., & Parrish, I. J. 2009, ApJ, 

704,211 
Bonafede, A., Feretti, L., Murgia, M., Govoni, F., Giovannini, G., Dallacasa, 

D., Dolag, K., & Taylor, G. B. 2010, A&A, 513, A30 
Brentjens, M. A. 201 1, A&A, 526, A9 
Brentjens, M. A., & de Bruyn, A. G. 2005, A&A, 441, 1217 
Carilli, C. L., & Rawlings, S. 2004, New Astronomy Review, 48, 979 
Carilli, C. L., & Taylor, G. B. 2002, ARA&A, 40, 319 
Churazov, E., Forman, W., Jones, C, Bohringer, H. 2003, ApJ, 590, 225 



Dewdney, R E., Hall, R J., Schilizzi, R. T., & Lazio, T. J. L. W. 2009, IEEE 

Proceedings, 97, 1482 
EnBlin, T A., & Vogt, C. 2003, A&A, 412, 373 
Feretti, L., & Giovannini, G. 2008, A Pan-Chromatic View of Clusters of 

Galaxies and the Large-Scale Structure, 740, 143 
Govoni, F, & Feretti, L. 2004, International Journal of Modern Physics D, 

13, 1549 
Guidetti, D., Murgia, M., Govoni, F, Parma, P., Gregorini, L., de Ruiter, 

H. R., Cameron, R. A., & Fanti, R. 2008, A&A, 483, 699 
Jarvis, M. J., & Rawlings, S. 2004, New Astron. Rev., 48, 1173 
Juett, A. M., Davis, D. S., & Mushotzky, R. 2010, ApJ, 709, L103 
Kelly, B., & McKay, T. 2005, AJ, 129, 1287 



Probing MHD Instabilities with Faraday Rotation 



11 



Kuijken, K. 2008, A&A, 482, 1053 

Krause, M., Alexander, P., Bolton, R., Geisbilsch, J., Green, D. A., & Riley, 

J. 2009, MNRAS, 400, 646 
Laing, R. A., Bridle, A. H., Parma, P, & Murgia, M. 2008, MNRAS, 391, 

521 
Massey, R., & Refregier, A. 2005, MNRAS, 363, 197 
Massey, R., Rowe, B., Refregier, A., Bacon, D. J., & Berg'e, J. 2007, 

MNRAS, 380, 229 
Million, E. T., Werner, N., Simionescu, A., Allen, S. W., Nulsen, R E. J., 

Fabian, A. C, Bohringer, H., & Sanders, J. S. 2010, MNRAS, 407, 2046 
Murgia, M., Govoni, E, Eeretti, L., Giovannini, G., Dallacasa, D., Fanti, R., 

Taylor, G. B., & Dolag, K. 2004, A&A, 424, 429 
Nagai, D., Vikhlinin, A., & Kravtsov, A. V. 2007, ApJ, 655, 98 
Padovani, R 2010, MNRAS, 1796 

Parrish, I. J., Quataert, E., & Sharma, P 2010, ApJ, 712, L194 
Parrish, I. J., Quataert, E., & Sharma, P 2009, ApJ, 703, 96 
Parrish, I. J., Stone, J. M., & Lemaster, N. 2008, ApJ, 688, 905 
Parrish, I. J., & Quataert, E. 2008, ApJ, 677, L9 
Parrish, I. J., & Stone, J. M. 2007, ApJ, 664, 135 
Parrish, I. J., & Stone, J. M. 2005, ApJ, 633, 334 
Pfrommer, C, & Jonathan Dursi, L. 2010, Nature Physics, 6, 520 



Quataert, E. 2008, ApJ, 673, 758 

Roettiger, K., Stone, J. M., & Burns, J. O. 1999, ApJ, 518, 594 

Ruszkowski, M., & Oh, S. R 2010, ApJ, 713, 1332 

Schilizzi, R.T., Alexander, P., Cordes, J.M., Dewdney, P.E., Ekers, R.D., 

Faulkner, A.J., Gaensler, B.M., Hall, P.J., Jonas, J.L., Kellermann, K.l. 

2007, SKA Memo 100, http://www.skatelescope.org 
Simard-Normandin, M., Kronberg, R P, & Button, S. 1981, ApJS, 45, 97 
Snowden, S. L., Mushotzky, R. F, Kuntz, K. D., & Davis, D. S. 2008, A&A, 

478,615 
Stepanov, R., Arshakian, T. G., Beck, R., Frick, P, & Krause, M. 2008, A&A, 

480, 45 
Strauss, M. A., Huchra, J. P, Davis, M., Yahil, A., Fisher, K. B., & Tonry, J. 

1992, ApJS, 83, 29 
Taylor, A. R. 2008, lAU Symposium, 248, 164 
Victor, J. D., Mechler, F, Ohiorhenuan, 1., Schmid, A., & Purpura K. P 2009, 

J. Neurophysiol. 102, 3414 
Wilman, R. J., et al. 2008, MNRAS, 388, 1335 
Young, C. A., Gallagher, P T., Ireland, J., & McAteer, R. 2005, American 

Geophysical Union SMSP 1 1, 07 



